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ABSTRACT 

A set of subroutines, written in 360 FORTRAN IV (level G or H), easy for the 
FORTRAN programmer to use, provides printed plots as part of normal output. These 
subroutines are simplified so that no choices rmist be made and generalized so that 
choices may be made by the programmer of the plotting characters, the scales, the ap- 
pearance of the grid, and other options . The FORTRAN routines generate ordinary out- 
put records (up to 132 characters) suitable for on- or off-line printing. 
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A SET OF 360 FORTRAN (LEVEL G OR H) SUBROUTINES FOR 
GENERATING PRINTED PLOTS 
by Lois T. Dellner 
Lewis Research Center 

SUMMARY 

A set of subroutines, written in 360 FORTRAN IV (level G or H), easy for the 
FORTRAN programmer to use, provides printed plots as part of normal output. These 
subroutines are simplified so that no choices must be made and generalized so that 
choices may be made by the programmer of the plotting characters, the scales, the ap- 
pearance of the grid, and other options. The FORTRAN routines generate ordinary out- 
put records (up to 132 characters) suitable for on- or off-line printing. 


INTRODUCTION 

Various versions (refs. 1 to 4) of this plotting system have been in use at Lewis 
Research Center since April 1962. One current version is a set of 360 FORTRAN sub- 
routines which generates ordinary FORTRAN output records (up to 132 characters) suit- 
able for printing on- or off-line. 

The system permits, but does not require, the programmer to choose the plotting 
characters, the scales, the grid-line spacing, etc. It handles single or multiple curves, 
prints true scales, and permits titles above or below the plot. The next section is a 
USERS MANUAL. 

The routines are presently in use at Lewis on an IBM 360-67 with no supervisory 
routine controlling printing. The programs are almost entirely machine independent, and 
the information provided in the section SYSTEM MANUAL is intended to make it very sim- 
ple for the systems programmer to make any changes required to adapt the plotting sys- 
tem to other machine configurations . The final section contains examples of plots obtain- 
able with these subroutines (see pp. 61 to 65). 



USERS MANUAL 


This system offers printed plots as part of normal output with a minimum of pro- 
gramming effort. The programmer writes CALL PLOTXY (for a single curve) or CALL 
PLOTMY (for multiple curves). The arguments, or call list, include the names of the 
arrays to be plotted and specify the number of points per curve and the number of curves. 
He may precede the call by writing a title to be printed above the plot. He may follow the 
call by writing a legend to be printed at the bottom of the plot. 

The plot or plots are printed as part of the regular output listing with no delay. No 
changes in his card-handling procedures nor special instructions for the operators are 
required. 

If he is using PLOTXY, the values of the variable to be plotted in the x-direction 
must be in sequence. If they are not, the subroutine SORTXY (which makes the neces- 
sary rearrangements) is supplied to be used before calling PLOTXY. For either PLOTXY 
or PLOTMY, if the size of the elements in (or the total range of) any array is not known 
to be within certain limits, the programmer calls the subroutine SHALE for each array 
before calling the plotting subroutine. SHALE will transform the array to suit PLOTXY 
and PLOTMY only if it is necessary. 

In addition to the minimum -effort use just described, the programmer may choose to 
use one or more of several options that permit him to control, for example, the appear- 
ance of the grid (by specifying the frequency of the grid-lines in either direction), the 
scale for either variable (by specifying the scale factor and a starting value) , the plotting 
character, etc. 

This section contains detailed instructions for the use of PLOTXY and PLOTMY and 
brief descriptions of SHALE and SORTXY. Some assistance in debugging users pro- 
grams completes this section. 

I. PLOTXY 

To get plotted output using PLOTXY, the corresponding pairs of ordinates to be 
plotted must be in two arrays. For an example, let us name the arrays XDOWN and 
YACROS and assume each is NPTS elements long. These names are chosen specifically 
to call the user’s attention to the fact that the x-direction is down the page. 

A. Call 

Write CALL PLOTXY (XDOWN, YACROS, KODE, P) 

XDOWN is the name of the array containing the values of the variable to be plotted 
on the x scale (down the page) . The elements in the array are restricted 
as follows: 
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(1) They must be in floating point. 


(2) The absolute value of each element must'*' be within permissible 

fi 

limits (approximately 10 < e < 10 ). 

2 

(3) They must be in order, either increasing or decreasing. 

YACROS is the name of the variable to be plotted on the y scale (across the page). 
The elements of the array are restricted exactly as in (1) and (2) above. 

KODE is the name of an integer . Many options are provided for the user of 

PLOTXY. Each has a number associated with it. The sum of the num- 
bers representing the options being used is KODE (0 < KODE < 119) . For 
your first plot, use KODE = 0. When KODE = 0, the starting-values and 
scale-factors in both directions are computed by PLOTXY. The other 
effects of KODE = 0 are shown in the example plot: 


(1) The plotting character is an asterisk. 

(2) The grid-line frequency is 10 x 10. 


(3) Nothing is printed to the left of the plot. 


P is an array. 

P(l) must contain NPTS (the number of points to be plotted) in floating 
point. 

If KODE = 0 or 64, there are no other requirements for P. 

For all other values of KODE, the requirements for the P array are dis- 
played (see section D. Using the Options, p. 6). 


*lf the element size of any array is unknown or out of range, write CALL SKALE 
(NPTS, A, KA) before calling PLOTXY (see III, p. 13). 

o 

If the array to be plotted on the x scale is not in order , write CALL SORTXY 
(XDOWN, YACROS, NPTS) before calling PLOTXY (see IV, p. 13). 
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B. Title 


The call for PLOTXY should be preceded by a WRITE statement causing a skip to 
a new page. This may be combined with a title to be printed above the plot. A repre- 
sentative example is: 

WRITE (6, 500) 

500 FORMAT (2H1 , 7 IX, 14HSAMPLE EXAMP LE/2H , 69X, 

17H JOLO SYSTEM PLOT) 

The result is shown on the plot on the facing page . 

C. Legend 

The call for PLOTXY should be followed by a WRITE statement causing a skip to a 
new page. This may be combined with a legend to be printed below the plot. An exam- 
ple is: 


WRITE (6, 502) 

502 FORMAT (72X, 16HF SUB 2 VS. THETA/72X, 
16H ( ) (RAD)/1H1) 

The result is shown on the plot on the facing page. 
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SAMPLE EXAMPLE 
JOLO SYSTEM PLOT 
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SUB2 VS. THETA 
— ) (RAD) 



D. Using the Options 

Each option has a number associated with it. The sum of the numbers repre- 
senting options chosen is KODE, the third argument in the call. 


TO CHOOSE : 

TO 

KODE 

ADD 

AND 

SUPPLY 

IN 

(FORMAT) 

KODE 0 
USES 

The plotting character . 

1 

Desired 

P(2) 

1H0 

* 

9 represents any acceptable FORTRAN 


charac- 




character except the minus sign. 


ter 

- 



The frequency of x grid-lines. 

2 

m 

P(3) 

Floating 

10 

They are printed every m line -spaces 




point 


below the first. If m = 0, only two 






will be printed, one above and one 






below the plot. 






The frequency of y grid-lines. 

4 

n 

P(4) 

Floating 

10 

They are printed every n positions 




point 


to the right of the first. If n = 0, 






only one is printed, at the left of 






the plot. 






The x scale. a 

16 

Three 

P(6), 

' 

Floating 

Scale 

The scaling parameters FX, repre- 


scaling 

P(7), 

point 

computed 

senting the star ting- value, and DX, 


param- 

P(8) 


by 

representing the scale -factor for one 


eters, 



plotting 

line space, must be whole numbers of 


KSX, FX, 



routine 

magnitude less than 10®. (0 < KSX < 6) 


and DX 




Seel.-E. Scaling. 






The y scale. a 



] 



The scaling parameters FY, repre- 

32 

Three 

P(9), 

Floating 

Scale 

senting the starting-value, and DY, 


scaling 

P(10), 

point 

computed 

representing the scale -factor for one 


param- 

P(ll) 


by 

print position, must be whole numbers 


eters, 



plotting 

of magnitude less than 10®. (0 < KSY < 6) 


KSY, FY, 



routine 

See I. -E. Scaling. 


andDY 




To print coordinates. 

64 

No re- 



No 

The coordinates of each point will be 


quire- 



print- 

printed on the line on which the point is 


ments 



out 

plotted. 







a To specify only FX, omit DX in P(8) and place 1. in P(5). 

To specify only FY, omit DY in P(ll) and place 2. in P(5). 

To specify FX and FY only , omit DX and DY and place 3. in P(5). 
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E. Scaling 


When the user wishes to specify his own scale in either direction, in addition to 
increasing KODE by 16 or 32, he must: 

(1) Choose his desired starting-value F, 

(2) Choose his desired scale-factor D - for one line-space if he is specifying the 

x scale, for one print position if he is specifying the y scale. 

(3) Determine a value of N such that: 

(a) F x 10 N is a whole number. 

(b) D x 10 N is a whole number. 

(c) 0 <N < 6. 

(4) Compute 6 - N. 

The integer calculated in step (4) is the first of three scaling parameters that must 
be available in the P array (in P(6) if scaling x, P(9) if scaling y) when the plotting sub- 
routine is called. F x 10^ is the second scaling parameter, and is placed in P(7) or 
P(10). F may be zero. D x 10^ is the third, and is placed in P(8) or P(ll). Note that 
D may never be zero, and when scaling x for PLOTMY, D must be positive. 

If the user wishes to supply only the starting value and leave the choice of scale- 
factor to the plotting system, he must do steps (1), (3) (a), (3)(e), and (4) above, and 

(a) Supply 1. in P(5) if option 16. 

(b) Supply 2. in P(5) if option 32. 

(c) Supply 3. in P(5) if both 16 and 32 are being used. 
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n. PLOT MY 


Although more than one curve can be plotted on the same grid with PLOTXY, only 
one plotting character will be used for all curves. PLOT MY provides a different plotting 
character for each curve. Options, similar to those in PLOTXY, are available but not 
required, except that option 1 must be used if there are more than six curves. However, 
a selection must be made from three Variations: DUPX - when more than one set of y 
values corresponds to the same set of x values; DUPY - when more than one set of x 
values corresponds to the same set of y values (this offers complete control of which 
variable is to be plotted in which direction); and NO DUP - when each set of x values has 
a corresponding set of y values. 

A. Call 

Write CALL PLOT MY (XDOWN, YACRQS, KKK, P) 

XDOWN is the name of the array containing the values of the variable to be 

plotted on the x scale (down the page) . The minimum DIMENSION of 
this array depends on the Variation selected. (See E . Variation Layout, 
pp. 10 to 12). The elements of this array are destroyed by PLOTMY. 
The elements are restricted as follows: 

(1) They must be in floating point, 

(2) They must be within permissible range. (Seel. PLOTXY, p. 2.) 

YACROS is the name of the array containing the values of the variable to be plot- 
ted on the y scale (across the page). The minimum DIMENSION of this 
array depends on the Variation selected, (See E. Variation Layout.) 

The elements of this array are restricted exactly as in (1) and (2) above. 

KKK is the name of an array. The first element must be KODE (as in 

PLOTXY, this is the sum of the option numbers), and the second ele- 
ment must be KN (the number of curves). The third element must be: 

(1) The number of points in one curve for DUPX and DUPY. 

(2) The number of points in the first curve for NO DUP. 

The remaining odd-numbered elements are only required for NO DUP, 
and the remaining even-numbered elements are only required for op- 
tion 1. (See E. Variation Layout.) The minimum DIMENSION of KKK 
must be 2 * KN + 2 or 14, whichever is greater . 


8 



P is the name of an array. The first element must specify the Variation 

selected: 

P(l) = 1. for DUPX 
P(l) = 3. for DUPY 
P(l) = 5. for NO DUP 

The remaining elements and the minimum DIMENSION of P are func- 
tions only of the options being used. (See E. Variation Layout. ) 

B. Title 

A page skip should occur before calling PLOTMY, exactly as for PLOTXY 
(p. 4). 

C. Legend 

A page skip should follow the CALL to PLOTMY, exactly as for PLOTXY 
(p. 4). 

D. Using the Options 


The use of the options is exactly the same as for PLOTXY (see I. -D. Using the 
Options, p. 6) with the following exceptions: 

(1) If option 16 is used, DX must be positive. 

(2) When option 64 is used, only the ordinate in the x-direction will be printed. 

(3) Unless option 1 is chosen, the plotting characters *, +, 0, X, =, 0 will be 
used for the first six curves. To use others, use option 1 and supply the desired 
plotting character for the first curve in KKK(4) , for the second curve in 
KKK(8), .... If more than six curves are to be plotted, option 1 must be used, 
and afl plotting characters supplied. To put a plotting character, for example, $, 
in KKK(4) , write 


DATA KKK(4)/1H$/ 
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E. Variation Layout 


DU PX 

[ ] means supply item in brackets if option N is being used. 



DIMENSIONS 

X(NPTS) 

Y(KN * NPTS) 

KKK •- (2 * KN + 2) or (14), whichever is greater 

(P«ll) , depending on use of options 2, 4, 

16, and 32. 


DUPY 



DIMENSIONS 

X(KN * NPTS) 

Y(KN * NPTS) 

KKK = (2 * KN + 2) or (14), whichever is greater. 

P(<11), depending on use of options 2, 4, 

16, and 32. 



NO DUP 


[ ]^. means supply item in brackets if option N is being used. 



DIMENSIONS 

X(NP1 + NP2 + NP3 + . . .) 

Y(NP1 + NP2 +NP3 + . . .) 

KKK « (2 * KN + 2) or (14), whichever is greater 

P(<11), depending on use of options 2, 4, 

16, and 32. 


HI. SKALE 


CALL SKALE (NPTS, A, KRSTR) 

This subroutine finds the largest absolute value of the NPTS elements of A and 
computes the characteristic of its logarithm to the base 10. If the characteristic 
K is -2 < K < 4, KRSTR is set to zero and control returns to the calling program. If 
K > 4 or < -2, each element of A is multiplied by a power (KRSTR) of 10 to transform 
the array to suit PLOTXY and PLOTMY. KRSTR is returned to enable the user to 
ReSToRe the array or to record how it has been altered. 

IV. SORTXY 


CALL SORTXY (V, W, NPTS) 

This subroutine rearranges the NPTS elements of the V array in order of increasing 
size. The elements of the W array are moved to maintain the original pair -relation; that 
is, if the fifth element of the V array is moved to the first position of V, the fifth element 
of W is moved to the first position of W. 

V. ERROR MESSAGES 

A. The message BAD LABELS printed below and to the left of the plot indicates that 
one or more of the x or y grid-labels may be incorrect. This can be caused by DX 
or DY too large or too small or by a label requiring more than the allotted nine print 
positions. 

B. When the system computes the scale factor for either array and finds the range is 
zero, calculations are stopped. For either PLOTXY or PLOTMY the plot terminates 
with a message: 

(1) For PLOTXY - N.G. followed by the contents of X(l), Y(l), X(2), Y(2), KODE, 
and P(l) 

or 

(2) For PLOTMY - N.G. followed by the contents of X(l), Y(l), X(2), Y(2), K(l), 
K(2) , and K(3). 
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C. If the values in the X array are not in order when PLOTXY is called, or if you use 
option 16 and a value of X lies outside your specified starting value, the message 

X OUT OF ORDER. I=NNNNN 

will be printed followed by the message described in (2) on the previous page, and the 
plot is terminated. 

D. In PLOTMY if the plotting character search fails, the message 

ERROR in K ARRAY 

is printed. This usually means DIMENSION of KKK is incorrect. 


SYSTEM MANUAL 

The objective of this section is to simplify the work of the systems programmer who 
implements and may have to modify this plotting system. To such a programmer the sys- 
tem consists of a set of FORTRAN subroutines (PLOTXY, PLOTMY, PISTUG, SORTXY, 
SKALE, and for level G only, a function subprogram AND) which are executed in a ma- 
chine with no supervisory printing control. A brief description and a complete listing 
are included for all routines . For PLOTXY, PLOTMY, and PISTUG, dictionaries of the 
FORTRAN variable names and comprehensive block diagrams are supplied. These are 
intended as systems debugging aids. Each dictionary is designed for use primarily as a 
cross-reference with the corresponding FORTRAN listing and block diagram. Within 
each dictionary, each variable name is identified in terms of its contribution to the pro- 
gram. Throughout these sections, the subroutine name is printed at the bottom of each 
page for assistance in rapid cross-referencing. 

The final discussion in this section is of the provisions made for changes that may be 
required if the plotting system is used with supervisory printing control. 

I. PLOTXY 

A. Description 

To plot one curve, the programmer will usually use PLOTXY. His call statement 
lists the names of the arrays to be plotted, the number KODE (an indicator of which op- 
tions are being used), and finally, an array whose first element is the number of points 
to be plotted and whose other elements are any additional data required by the choice of 
option. 

PLOTXY writes a blank line (FORMAT 2H ) to initiate the plot, and calls the 
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auxiliary subroutine PISTUG (p. 45) to compute scaling parameters for either or both 
scales and to check label sizes . Communication with PISTUG is through 

COMMON/JO LO/F, DX, TIINX,N, LABOUT, KPWR, KFD 

COMMON/JO LO/XYX, FORY, STUG, TONLY, KSW64 

Records are then written for the plot beginning with a line of 11 y grid-labels. For 
each line of the plot any required horizontal and vertical grid-line characters as well 
as correctly positioned plotting characters for all points in that line are written. An 
x grid-label is printed every tenth line . When all points have been plotted, the plot is 
terminated (with an x grid-line at the next x grid-label, followed by a line of y-labels 
and a final blank line (FORMAT 2H )) . 

Since the values in the X array are required to be in order (up or down) , no sorting 
is required and the user’s arrays are undisturbed. 


PLOTXY 
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B. Program Listing 


SUBROUTINE PLOTXY ( X , Y ,K , P ) 

C FIV-DCS TO 360 (G OR H) 

LOGICAL X YX , FURY , STUG , TONLY , XGL » LS , K SW64 

DIMENSION FYLAB ( 6 ) , YLABEL(ll), A(104) , FXLABSI6), FXL ABM (7 ) 
DIMENSION X ( 1 ) ,Y( 1) ,P( 1) 

COMMON/JOLO/F,DX,TLINX,N,LABOUT, KPWR,KFD 
COMMON/ JOLO/XYX, FORY, STUG, T0NLY, KSW64 
EQUIVALENCE ( FYL AB ( 5 ) , IYLAB) , ( FXL AB , I FXL AB ) 

DATA MASK 1 , MASK2 , MASK4, MASKS , MASK 16 , MA SK 32 , MASK 64 / 

1 1,2,4,8,16,32,64 / 

DATA FYLABt 1 ) , FYL AB ( 2 ) , FYLAB ( 3 ) , FYLAB ( 4 ) , F YL AB5 , FYL AB ( 6 ) / 

1 4H ( 2 H , 4H ,20, 4HX , , 4H11F , 4H10.0, 4H ) / 

DATA BLANK , XGR I D , Y GR I D /1H ,1H-,1H1 / 

DATA PCSTD/1H*/ 

DATA RMARK/1H= / 

DATA FXLA3S ( 1 ) ,FXLA8S(2) , FXLA8SI 3) , FXLAB I , F XL ABS ( 5 ) , FXL A8S ( 6 ) / 

1 4H(2H , 4H ,18, 4HX , , 4HF9.0, 4H,104, 4HA1 ) / 

DATA FXLABMI 1 ) , FXL ABM ( 2 ) , FXL ABM( 3 ) , FXLABM ( 4 ) , F XL ABM ( 6 ) , FXL ABM ( 7 ) / 
1 4H ( 2 H , 4H ,4X, 4H,2F7, 4H.3, , 4H,104, 4HA1 ) / 

C 

100 WRITE (6,500) 

102 KODE=K 
N=P ( 1 ) 

L AB0UT=1 
FXLAB = FXLABI 
LS = .FALSE. 

F Y L A B ( 5 ) = F Y L A B 5 
KSW64 = .FALSE. 

110 PC=PCSTD 

112 I F ( ( AND ( KOBE , MASK 1 ) ) . GT, 0 . ) PC = P(2) 

114 M=10 

116 I F ( ( AND ( KODE , MASK 2 ) ) . GT. 0 . ) M=P ( 3 ) 

117 IF ( M.EO.O )M=1000 

118 NY-10 

120 IF( (AND(K0DE,MASK4) ) .GT.O. )NY = P(4) 

121 IF(NY.EQ.O) NY= 1 000 

122 I F ( ( AND ( KODE , MASK64 ) ) .GT .0. )KSW64=.TRUE. 

124 I F ( (AND (KODE, MASKS) ) .GT.O. ) GO TO 1240 
126 GO TO 140 

C 

1240 WRITE (6,530) 

C 

140 XYX=. FALSE. 

142 FORY= . TRUE . 

144 STUG*. FALSE, 

146 TONLY=. FALSE. 
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148 I F ( ( AND(K0DE,MASK32 ) .LE.O. ) ) GO TO 172 

151 ST!JG=.TRUE. 

152 KSY=P ( 9 ) 

154 PWR10Y=10.**(KSY-6) 

156 FY -P( 10)*PWR10Y 
158 F = FY 
C 

160 IF(P(5).GE.2.) GO TO 172 
162 TONL Y= . TRUE . 

164 DY= P ( 1 1 ) *PWR 10 Y 
166 DX= DY 

172 CALL PISTUG(Y) 

173 IF( DX. EO. 0. ) GO TO 700 

174 FY = F 
176 DY=DX 

180 I F ( KSW64 ) KPWRY = KPWR 
190 I YLAB=I YLAB+KFD 
C 

200 XYX = .TRUE . 

202 FORY=. FALSE. 

204 STUG=. FALSE. 

206 TONLY=. FALSE. 

208 TLINX=55*< l+N/35) 

C 

210 IF( ( AND(K0DE,MASK16) .LE.O. )) GO TO 232 

213 STIJG= .TRUE . 

214 KSX = P ( 6 ) 

216 PWR 10X=01 0 . *4 ( KSX-6 ) 

218 F X= P ( 7 ) *PWR 1 OX 
220 F=FX 
C 

222 IF(MOO( I F I X ( P ( 5) ) ,2) .EQ.l) GO TO 232 
224 TONL Y=. TRUE. 

226 DX =P ( 8 ) *PWR10X 
232 CALL PISTUG(X) 

IFIDX.EO.O. ) GO TO 700 
234 FX=F 

240 I F ( KSW64 ) KPWRX = KPWR 
242 IFXLAB = IFXLAB +KFD 
FXL ABS ( 4 ) = FXLAB 
FXLABM ( 5 ) = FXLAB 
C 

250 I F ( . NOT . K S W64 ) GO TO 264 

252 KOUT X = -KPWR X 

254 KOUTY=— KPWRY 

256 F10X=10.**KPWRX 

258 F10Y=10.*4KPWRY 

260 WRITE (6,502) KOUTX,KOUTY 


PLOTXY 
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264 DO 278 1=1,11 

266 TEMP = FY+ FLOAT ( I- — 1 ) *DY* 10 • 

268 A T E M P = ABS(TEMP) 

270 IF (ATEMP.LT.l.E-7) TEMP = 0. 

272 IF ( ATEMP .GE . 1 . E+7 ) LAB0UT=2 
278 Y LABEL ( I )=TEMP 
300 KSYLAB =1 

302 WRITE ( 6 , FYLAB ) ( YL ABE L ( I ) , I = 1 , 1 1 ) 
304 GO TO (306,700) , KSYLAB 
C 

306 KSYLAB =2 
310 LC TR = 0 
NCTR=1 
KGUT=1 
KQU I T = 1 
C 

320 I F ( MOD { LCTR , M ) >328,322,328 
322 XGL = . TRUE • 

AFI LL= XGRID 
GO TO 330 
328 XGL =. FALSE. 

AF ILL-BLANK 
330 DO 332 1=2,104 
332 A (I) = AFILL 
334 DO 336 1=2, 104, NY 

336 A ( I ) = YGRID 
A(l) =BL ANK 

338 GO TO (340,400) ,KOUT 
C 

340 KX =( X(NCTR)-FX )/DX +*5' 

342 I F (KX-LCTR >630,350,600 

350 KY= ( Y(NCTR)-FY)/DY+.5 

351 LS= .TRUE • 

352 TPC = PC 

353 KYL = KY+2 

354 IF(KY.LT.O) GO TO 360 
356 I F ( K Y .G T • 1 0 1 ) GO TO 364 
358 GO TO 370 

360 KYL= 1 
362 GO TO 366 
364 K YL= 1 04 
366 TPC=RMARK 
C 

370 A (KYL) =TPC 
380 I F ( .NOT .KSW6 4 > GO TO 386 
382 XOUT = X ( NC TR ) / F 1 0 X 
YOUT = Y ( NCTR ) / F 10Y 
C 

386 I F ( NCT R . GE • N ) GO TO 392 
388 NCTR = NC TR+ 1 
390 GO TO 340 
C 

392 KOUT = 2 
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400 I F ( MOD ( LCTR , 1 0 ) .GT.5) M = 10 

401 I F { XGL.AND. (MOD (LCTR* 10).BQ.0) )K0UIT=2 
C 

600 JX - 1 

602 I F ( MOD ( LCTR , 10 ) *NE • 0 ) GO TO 613 
604 JX = JX+2 

606 XLABEL =FX+ FLOAT (LCTR)*DX 
608 TEMP -ABS< XLABEL) 

610 I F ( TEMP .GE • 1 • E+7 ) LAB OUT = 2 

612 I F ( TEMP . LT. 1 .E-7 ) XLABEL-O. 

613 I F ( .NOT .KSW64 ) GO TO 6240 
C 

620 3 F ( .NOT *LS ) GO TO 6240 

621 JX = JX+1 

622 LS = .FALSE. 

C 

6240 GO TO (6250,6260,6270,6280) , JX 
6250 WRITE <6, 504) (A (I) , 1=1 ,104 ) 

GO TO 614 
C 

c 

6260 WRITE (6, 505) XOUT , YOUT , ( A ( I) , I = 1,1.04) 

GO TO 614 
C 

6270 WRITE (6 t FXLABS) XLABEL, (A(I ) , I = 1,104) 

GO TO 614 
C 

6280 WR I TE ( 6 , F XL ABM ) XOUT , YGUT , XL ABEL , ( A ( I) , I = 1,104) 
C 

614 LCTR = tCTR+1 

616 GO TO ( 320,302 ) , KQU I T 
C 

630 LABOUT = 4 

700 GO TO (710, 702, 706, 704), LABOUT 
702 WRITE (6,506) 

GO TO 720 


C 


704 

706 

708 

710 

720 

500 

502 

504 

505 

506 
508 
510 
520 
530 


WRITE ( 6,520) NCTR 

WRITE (6,508) < X .( I ) , Y ( I ) ,1=1,2) , K , P ( 1 > 


GO TO 720 
WRITE <6 ,510) 
RETURN 
FORMAT (2.H ) 

FORMAT (2H , 

FORMAT (2H , 

FORMAT <2H , 

FORMAT (2H , 

FORMAT (2H , 

FORMAT (2H ) 

FORMAT ( 2H , 
FORMAT ( 2H 
END 


6X , 3HX*E, 12, IX, 4H Y*E ,I2 ) 

2 7X , 104 A 1 ) 

5.X , F6.3, IX, F6.3, 9X,104A1 ) 

3X , 1 OH BAD LABELS ) 

5H N.G • , 4G20. 8 , I 6 ,F8 « 2 ) 

18HX OUT OF ORDER. I = , 15) 

, T 100 , 1 OPT ION 8 NO LONGER AVAILABLE* 
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C. Block Diagram 
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D. PLOTXY Dictionary 


A FILL 

ATEMP 

BLANK 

DX 

DY 

F 

FI OX 

F10Y 

FORY 

FX 


Array in which each line of the plot proper is constructed. It is 
initialized to 104 blanks or 104 minuses (the horizontal grid-line char- 
acter). The vertical grid-line characters (l's) are placed in position 2 

i-U 

and every NY position thereafter. Finally, the plotting characters 
(if any) are positioned. 

Temporary storage for the character ( or -) with which the A array 
is initialized. 

Temporary storage for the absolute value of each y grid-label. 

Internal representation for a blank printing character. It is used to 
initialize the A array for lines that are not x grid-lines. 

Vertical scale factor (per printing line) . It is either calculated by 
PISTUG based on the range and length of the X array or it is calculated 
from values supplied by the user of option 16 in P(6) and P(8) as 
DX = P(8) • 10 p ( 6 >- 6 . It is used to calculate KX and XLABEL. Also 
see DX, page 52. 

Horizontal scale factor (per print position) . It is either calculated by 
PISTUG based on the range and length of the Y array, or it is calcu- 
lated from values supplied by the user of option 32 in P(9) and P(ll) as 
DY = P(ll)‘ loW- 6 . It is used to calculate KY and the values of 
Y LABEL. Also see DX, page 52. 

See COMMON/JOLO/ - Dictionary, page 52. 

10 KPWRX . Divisor of all values in the X array before they are printed 
at the left of the plot when option 64 is used. 

KPWPV 

10 . Divisor of all values in the Y array before they are printed 

at the left of the plot when option 64 is used. 

See COMMON/JOLO/ - Dictionary, page 52. 

Starting value for the vertical scale. It is either calculated by PISTUG 
based on the range, length, and XI (see p. 51 ) of the X array, or it is 
calculated from values supplied by the user of option 16 in P(6) and 
P(7) as FX - P(7) • 10 P ^" 6 . It is used to calculate KX and XLABEL. 
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FXLAB 

FXLABI 

FXLABM 

FXLABS 

FY 

FYLAB 

FYLAB5 

I 

I FXLAB 


A word initialized to F9.0 and altered to F9.d by using KFD(see p. 52). 
FXLAB is then moved into the FORMAT arrays (FXLABS and FXLABM) 
used for lines that have an x-label. It is equivalenced to IFXLAB. 

Initial value for FXLAB. 

Variable FORMAT array used to print a line when JX = 4 (see p. 25) . 

It is of the form 


(2H , 4X, 2F7.3, F9.d, 104A1) 

where 0 < d < 6. The value of d is inserted by the program (using 
FXLAB) after the call to PISTUG to process the X array. See KFD 
(P- 52). 

Variable FORMAT array used to print a line when JX = 3 (see p. 25). 
It is of the form 


(2H , 18X, F9. d, 104A1) 

where 0 < d < 6 . A value of d is inserted by the program (using 
FXLAB) after the call to PISTUG to process the X array. See KFD 
(p. 52). 

Starting value for the horizontal scale. It is either calculated by 
PISTUG based on the range, length, and XI (see p. 51) of the Y array, 
or it is calculated from values supplied by the user of option 32 in P(9) 
and P(10) as FY = P(10) • 10 P ^" 6 . It is used to calculate KY and the 
values of Y LABEL. 

Variable FORMAT array with which the line of y grid-labels is writ- 
ten. It is of the form 


(2H , 20X, llFlO.d) 

where 0 < d < 6. A value of d is inserted by the program (as KFD) 
after the call PISTUG to process the Y array. See KFD (p. 52). 

Initial value for FYLAB(5) . 

Index of loop forming y grid-labels and of loop initializing print line in 
A array. 

Equivalent to FXLAB. 
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IYLAB 

JX 

K 

KFD 

EODE 


KOUT 

KOUTX 

KOUTY 

KPWR 

KPWRX 

KPWRY 

KQUIT 

KSW64 

KSX 

KSY 


Equivalent to FYLAB(5) . 

Initialized to 1, increased by 1 if option 64 is being used and a point 
has been placed in the line, increased by 2 if the line to be printed has 
an x-label, this switch value may be 1, 2, 3, or 4. It is used to con- 
trol the selection of the correct WRITE statement for each line . 

Third argument in the call list of this subroutine. It contains KOBE. 

See COMMON/JOLO/ - Dictionary, page 52. 

Sum of the option numbers. It is supplied by the programmer in the 
third argument of the call. Branching on single bit positions using 
MASK1, MASK2, MASK4, MASKS, MASK16, MASK32, and MASK64 
serves to identify the options being used. Option 8 has been discon- 
tinued but a message is printed if MASK8 causes a branch. 

Switch initialized to 1 and set to 2 after all points have been processed. 

Value in the special heading of option 64 for the x coordinates. It is 
the negative of KPWRX. 

Value in the special heading of option 64 for the y coordinates . It is 
the negative of KPWRY. 

See COMMON/JOLO/ - Dictionary, page 52. 

KPWR for the X array . 

KPWR for the Y array. 

Switch initialized to 1 . It is set to 2 just before the last line (a labeled 
x grid-line) of the plot is printed. 

See COMMON/JOLO/ - Dictionary, page 52. 

Fixed point form of the value supplied in P(6) by the programmer 
using option 16. P(7) and P(8) are multiplied by to get FX 

and DX. 

Fixed point form of the value supplied in P(9) by the programmer 
using option 32. P(10) and P(ll) are multiplied by to get FY 

and DY. 
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KSYLAB 

KX 

KY 

KYL 

LABOUT 

LCTR 

LS 

M 

MASK1 

MASK2 

MASK4 

MASK8 

MASK16 

MASK32 

MASK64 

N 

NOTE 


Switch initialized to 1 and set to 2 after the first line of y grid-labels 
has been written. 

Line number of a particular point. It is computed from 

[(x. - FX)/DX] rounded. The first x grid-line is line number zero. 

"Printing position” of a particular point. It is computed from 
[(y. - FY)/DY] rounded. The first y grid-line character in "printing 
position" zero is actually the second element of the A array. 

Value of KY + 2 if the y -coordinate falls on the plot; 1 if KY is nega- 
tive; 104 if KY > 101. 

See COMMON/JOLO/ - Dictionary, page 52. 

Line counter. It is set to zero for the first x grid-line and is stepped 
by 1 after each line is printed. 

Logical variable set to TRUE whenever a point has been found for the 
current line . If option 64 is in effect, LS is used to control a branch 
(to modify JX) and is set to FALSE immediately after the branch. 

Frequency of the x grid-lines. It is either Supplied in P(3) by the pro- 
grammer if option 2 is used or set equal to 10. If the programmer 
specifies zero, M is set to 1000 and only the first and last x grid-lines 
will be printed. 


>-See KODE, page 25. 


See COMMON/JOLO/ - Dictionary, page 53. 

Counter for the number of points that have been processed. It is ini- 
tialized to 1 and tested for equality to N. It is increased by 1 each 
time the test fails. 
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NY 

P 

PC 

PCSTD 

PWR10X 

PWRIOY 

RMARK 

STUG 

TEMP 

TLINX 

TONLY 

TPC 

X 

XGL 

XGRID 

XLABEL 

XOUT 

XYX 


Frequency of the y grid-lines . It is either set to 10 or supplied in P(4) 
by the programmer using option 4. If it is supplied as zero, it is re- 
placed by 1000 and only the left-most y grid-line will be printed. 

Array name, the fourth argument in the call list of this subroutine. 

P(l) contains a number of points to be plotted. Further contents are 
prescribed by the options being used. 

Plotting character. It is either supplied by the programmer using op- 
tion 1 in P(2) or by the program as PCSTD. 

Standard plotting character (*) . 

Value of 10 KSX “ 6 . 

Value of lO^- 6 . 

Out-of-range character (=). 

See COMMON/JOLO/ - Dictionary, page 53. 

Temporary storage. 

See COMMON/JOLO/ - Dictionary, page 54. 

See COMMON/JOLO/ - Dictionary, page 54 . 

Current plotting character . 

Array name, the first argument in the call list of this subroutine, con- 
taining the values of the variable to be plotted on the vertical scale 
(down the page) They must be in sequence, either increasing or de- 
creasing. 

Logical variable set to TRUE whenever the value of LCTR is evenly 
divisible by M. 

x- grid -line character, a minus sign. 

Current value of the x label. It is computed only for every tenth line 
from FX + LCTR - DX. 

Value of Xj modified by FI OX for option 64 printout. 

See COMMON/JOLO/ - Dictionary, page 54. 
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Y 


YGRID 

YLABEL 

YOUT 


Array name, the second argument in the call list of this subroutine, 
containing the values of the variable to be plotted on the horizontal 
scale (across the page) . 

y grid-line character - the digit 1. 

Array holding the 11 labels for the y axis. 

Value of Y. modified by F10Y for option 64 printout. 
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II. PLOT MY 


A. Description 

Although more than one curve can be plotted on the same grid using PLOTXY, they 
will all use the same plotting character. When a different plotting character for each 
curve is desired, the programmer uses PLOT MY. The CALL statement lists the 
names of the arrays to be plotted; an array containing the number KODE , the number 
of curves, and possibly other information - depending on the options used and the 
Variation chosen; and finally an array whose first element indicates the Variation 
chosen and whose other elements are any additional data required for the options being 
used. 

The records written by PLOTMY are identical to those written by PLOTXY 
(p. 14), and the auxiliary subroutine PISTUG is used in the same way. 

Because the X array is not required to be in order, a search for the next point to 
be plotted occurs for every point. The search technique used avoids the use of addi- 
tional storage by using the X array itself in such a fashion that the values in the array 
are destroyed. The Y, KKK, and P arrays are undisturbed. 
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B. Program Listing 


SUBROUTINE PLOT MY ( X , Y , K , P ) 

C DC 4 TO 360 

LOGICAL XYX , FORY f STUG, TQN'LY ,XGL , L S , KS W64 

DIMENSION FYLAB(6) ,YLABELlll) , AI104), XLAB(7) f XLAB6M7) 
DIMENSION KPCSTD (6 ) 

DIMENSION X( 1 > ,Y(1) ,PU ) ,K( 1,) 

COMMON/ JGLO/ F , DX ,TL I NX ,N,L ABOUT, KPWR,KFD 
COMMON/ JOLO/XYX, FOR Y , STUG, TONLY, KSW64 
EQU I VALENCE { X L A B 5 , I XL ABB) 

EQUIVALENCE ( FYLAB ( 5 ) , I YL AB ) 

EQUIVALENCE (KPC,TPC ) 

DATA BLANK, XGRTD,YGRID /1H ,1H-,1H1 / 
data RMARK/ 1H= / 

DATA KPCSTD /1H*,1H+,1H0,1HX,IH=,1HQ / 

DATA XL AB 5 I/4HF9 • 0/ 

DATA XL AB ( 1 ) , XLAB ( 2 ) , XLAB ( 3 ) , XL AB (4) , XL AB (6) ,XLAB(7)/4H (2H ,4H , 

I ,4H18X, ,4H , 4H, 1 04 , 4HA 1 ) / 

DATA XL AB64 ( 1 ) , XL AB64 ( 2 ) ,XLAB64(3> ,XLAB64(4) , XL AB64 ( 6 ) , XL AB6A- ( 7 ) / 
14H ( 2H ,4H , 12, 4HX, F, 4H6.3, ,4H, 104, 4HA1 ) / 

DATA MASK 1 , MASK 2 ,MASK4 , MASK8 , MASK 1 6 , MA SK 32 , MASK64 / 

1 1,2,4,8,16,32,64 / 

DATA FYLAB ( 1 ) , FYLAB ( 2 ) , FYLAB ( 3 ) , FYLAB ( 4 ) , F YL AB5 , FYL AB < 6 ) / 

1 4H(2H , 4H ,20, 4HX , , 4H11F , 4H10.0, 4H ) / 

C 

100 WRITE 16,500) 

102 K0DE=K(1) 

KN=K ( 2 ) 

NPTS=K ( 3 ) 

L ABOUT = 1 
LS = .FALSE. 

FYLAB (5)=FYLAB5 
KSW64 = .FALSE. 

XL AB5 = XL AB 5 I 
KTL-1 
K S W I = 1 

11000 I F ( P ( 1 ) -2 • 5 ) 11002,11002,11010 
11002 K$W 1=2 
11004 KTL=KN 
11006 GO TO 11034 

II 010 IF(P(1) -4. 0)110 2 0,11020, 11012 

11012 N P T S T = 0 

11014 DO 11016 1=1, KN 

11016 NPTST=NPT$T+K i 2*1+1 ) 

11017 TLINX=55*< 1+NPTST/ ( 35*KN ) ) 

11018 GO T0112 
11020 KTIMES=KN~1 

11022 DO 11032 I =1 , KT I MES 
11024 MM=I*NPTS 
11026 K ( 2*1+3 )=NPTS 
11028 DO 11032 11=1, NPTS 
11030 L = MM+ 1 1 
11032 Y ( L ) = Y ( I I ) 
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11034 NPTST=KN*NPTS 
11036 TL I NX -55* ( l + NPTS/35 ) 

112 I F ( { AND ( KODE ? MASK 1 ) ) . GT • 0 . ) GO TO 114 

11112 DO 11113 1=1 t KN 

11113 M2*I+2)=KPC$TD( I ) 

114 M= 10 

116 IF ( ( AMD (KODE .,KASK2.) ).GT.O. )K=P( 3) 

117 IF ( M * E 0 • G ) M = I 00 0 

118 N Y = 1 0 

120 IF ( « AND ( K ODE » M A S K4 ) ) * GT • 0 • 1 NY = P(4> 

121 IF(NY.EQ.O) NY=1 000 

122 I F i ( AND ( KODE , MASK64 ) ) .GT .0 . )KSW64=. TRUE . 

124 I F ( ( AND ( KODE * MASKS ) ) • GT. 0. ) GO TO 1240 

125 GO TO 140 
C 

1240 WRITE ( 6 1 530 ) 

C 

140 XYX=. FALSE. 

142 FORY = .TRUE . 

144 STUG=. FALSE. 

146 TONLY=. FALSE. 

C 

148 I F ( ( AN D ( KODE f M A SK 32 ) . L E • 0 * ) ) GO TO 172 

151 $TUG = .TRUE. 

152 KS Y=P ( 9 ) 

154 PWR10Y=10.** (KSY-6) 

156 FY =P( 10>*PWR10Y 
158 F = FY 
C 

160 IF ( P ( 5 ) .GE .2 • ) GO TO 172 
162 TONL Y= * TRUE . 

164 DY= P( ID^PWRIOY 
166 DX= DY 
C 

172 N-NPTST 
11172 CALL PISTUG(Y) 

173 IF( DX.EO.O . ) GO TO 700 

174 FY=F 
176 DY=DX 

190 IYLAB=IYLAB+KFD 
C 

200 XYX=. FALSE. 

202 FORY=. FALSE. 

204 ST'UG=. FALSE. 

206 TONLY= .FALSE • 

C 
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210 I F ( { AND(K0DE,MASK16) .LE.O. ) ) GO TO 232 

213 STUG= . TRUE • 

214 KSX = P ( 6 ) 

216 PWR10X=10.**(KSX-6) 

218 F X = P(7)*PWR10X 
220 F = FX 
C 

222 I F t MOD ( I F I X ( P ( 5 ) ) ,2 ) • EQ « 1 ) GO TO 232 
224 TONL Y= ♦ TRUE • 

226 DX =P(8)*PWR10X 

232 IFCKSWI • EQ • 2 ) N=NPTS 
11233 I L I M=N 

233 CALL PISTUG(X) 

I F ( DX • EQ .0 • ) GO TO 700 

234 FX=F 

I XLAB5=I XL AB5+KFD 
X L A 8 ( 5 ) = XLAB5 
XLAB64( 5 ) = XLAB5 
C 

250 IF( .N0T.KSW64JG0 TO 262 
252 KOUTX=-KPWR 
256 F10X=10.**KPWR 
260 WRITE ( 6 , 502 ) KOUTX 
C 

262 TDY = DY* 1 0 • 

264 DO 278 1=1,11 

266 TEMP = FY+ FLOAT ( I~1 )*TDY 

268 ATEM P= ABS (TEMP) 

270 IF (ATEMP.LT.l.E-7) TEMP = 0. 

272 IF ( AT EMP . GE • 1 • E + 7 ) LA60UT=2 
278 YLABEL ( I ) =TEMP 
300 KSYLAB =1 

302 WRITE ( 6 , F YL AB ) ( YLABEL (I)M =1 r 11 ) 

304 GO TO (306,700) , KSYLAB 
C 

306 KSYLAB =2 
310 LC TR=0 
NCTR=1 
KQUT = 1 
KQU I T = 1 
C 

320 I F (MOD ( LCTR ,M ) >328,322,328 
322 XGL = • TRUE • 

324 A FILL- XGRID 
GO TO 330 
328 XGL = . FALSE. 

AFILL=BLANK 
330 DO 332 1=2,104 
332 All) = AFILL 
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334 DO 336 1=2 ,104, NY 

336 A ( I ) = YGRID 
A { 1 ) =BL ANK 

33B GO TO (340*400) *KOUT 
C 

340 XMIN=1.E15 
11340 I M I N=1 

11342 DO 11350 1 = 1 1 H IM 

11344 IF ( XMIN— X ( I ) ) 11350* 1 1350* 1 1346 

11346 XM I N-X (I) 

11348 I M I N= I 
11350 CONTINUE 

341 KX =(Xt IMIN)-FX)/DX +.5 

342 I F { K X-LCTR ) 630 * 350 * 600 
350 LS=.TRUE. 

11352 X < I MI N ) = 1 # E l 5 

380 I F ( *NOT #KS W64 ) GO TO 11400 
382 XOUT = X M IN / F 1 0 X 
11400 DO 370 I M= 1 *KTL 
11402 LL= I.MIN+ ( I M -1)*NPTS 
11404 K Y= ( Y ( LL ) -F Y ) / DY+ *5 
11420 I F ( KS WI • EO • 2 ) GO TO 11440 
11422 IK=0 

KLAST=2^KN+1 

11424 DO 11430 IL=3*KLAST*2 
11426 I K = I K + K ( I L ) 

11428 I F C IK — IMIN) 11430*11436*11436 

11430 CONTINUE 

11432 L ABOUT = 5 

11434 GO TO 700 

11436 KPC=K( IL+1 ) 

11438 GO TO 353 
11440 KPC=K (2*1 M+2 ) 

3 53 KYL = KY4-2 

354 IF(KY.LT.O) GO TO 360 

356 IF ( K Y*GT • 101 ) GO TO 364 

358 GO TO 370 

360 KYL= 1 

362 GO TO 366 

364 K YL= 104 

366 TPC=RMARK 

370 A (KYL) =T PC 

386 IF(NCTR.GE.ILIM) GO TO 392 
388 NC TR=MCTR+ 1 
390 GO TO 340 
C 

392 KOUT = 2 
C 

400 I F ( ( MOD ( LCTR * 10 ) ) .GT.5)M = 10 

401 IF( XGL.A NO . ( MOD ( LC TR , 1 0 ) . E 0 . 0 ) ) KOU IT =2 
C 
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600 JX = 1 

602 I F ( MOD ( LCTR * 1 0 ) • NE • 0 ) GO TO 613 
604 JX = JX + 2 

606 XLABEL =FX+FL OAT ( LCTR ) *DX 
608 TEMP = ABS ( XL ABEL ) 

610 I F ( TEMP.GE.l.t+7)LAB0UT=2 

612 I F ( TEMP. LT • 1 • E-7 ) XL ABEL =0* 

613 I F ( .N0T.KSW64) GO TO 6240 
C 

620 IF(.NOT.LS) GO TO 6240 

621 JX = JX+1 

622 LS = .FALSE. 

C 

6240 GO TO ( 6250 * 6260 * 62 70 *62 80 ) * JX 
6250 WRITE (6*504) ( A ( I ) * I =1 * 104 ) 

GO TO 614 
C 

6260 WRITE (6*505) XOUT * ( A (I) * I = 1*104) 

GO TO 614 
C 

6270 WR I T E ( 6 * XL AB ) XL AB EL * ( A ( I ) * I = 1,104) 

GO TO 614 
C 

6280 WRITE (6*XLAB64) XOUT * XL ABEL * ( A ( I) ,1 = 1*104) 

C 

614 LCTR = LCTR+1 

616 GO TO (320*302 ) *KOUIT 
C 

630 LABOUT= 4 

700 GO TO (710*702*706,704*712) ,LABOUT 
702 WRITE (6*506) 

GO TO 720 
C 

704 WRITE(6*520) NCTR 

706 WR ITE (6*508) ( X ( I ) * Y ( I ) * I = 1 * 2 ) * ( K ( J ) * J= 1 * 3 ) * P ( 1 ) 
708 GO TO 720 
712 WRITE(6*512) 

710 W R I T E (6*510) 

720 RETURN 

500 FORMAT ( 2H ) 

502 FORMAT ( 2H * 1 3X , 3HX*E * I 2 ) 

504 FORMAT ( 2H *27X* 104A1 ) 

505 FORMAT ( 2H * 1 2 X * F6 . 3 * 9 X , 104 A 1 ) 

506 FORMAT ( 2H *3X*10HBAD LABELS ) 

508 FORMAT ( 2H *5H N . G . * 4G20 . 8 * 3 I 7 * F 8 . 2 ) 

510 FORMAT ( 2H ) 

512 FORMAT ( 2H * 1 6HERR0R IN K ARRAY ) 

520 FORMAT ( 2H *18HX OUT OF ORDER. I=*I5) 

530 FORMAT ( 2H * T 1 00 * • OPT I ON 8 NO LONGER AVAILABLE* ) 
END 
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C . Block Diagram 
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D. PLOT MY Dictionary 


A FILL 

ATEMP 

BLANK 

DX 

DY 

F 

F10X 

FORY 

FX 

FY 


Array in which each line of the plot proper is constructed. It is ini- 
tialized to 104 blanks or 104 minuses (the horizontal grid-line charac- 
ter). The vertical grid-line characters (l’s) are then placed in posi- 

j-L 

tion 2 and every NY position thereafter. Finally, the plotting char- 
acters (if any) are positioned. 

Temporary storage for the character ( or -) with which the A array is 
initialized. 

Temporary storage for the absolute value of each y grid-label. 

Internal representation for a blank printing character. It is used to 
initialize the A array for lines that are not x grid-lines . 

Vertical scale factor (per printing line) . It is either calculated by 
PISTUG based on the range and length of the X array or it is calculated 
from values supplied by the user of option 16 in P(6) and P(8) as 
DX = P(8) • 10 P ^~®. It must be positive. Also see DX, page 52. 

Horizontal scale factor (per print position) . It is either calculated by 
PISTUG based on the range and length of the Y array, or it is calculated 
from values supplied by the user of option 32 in P(9) and P(ll) as 
DY = P(U).10 P < 9 H Also see DX, page 52. 

See COMMON/JOLO/ - Dictionary, page 52. 

KPWRX 

10 . Divisor of all values in the X array before they are printed 

on the left of the plot when option 64 is used. 


See COMMON/JOLO/ - Dictionary, page 52. 

Starting value for the vertical scale. It is either calculated by PISTUG 
based on the range, length, and XI (see p. 51) of the X array, or it is 
calculated from values supplied by the user of option 16 in P(6) and P(7) 
as FX = P(7) • io P ^“®. It is used to calculate KX and XLABEL. 

Starting value for the horizontal scale. It is either calculated by 
PISTUG based on the range, length, and XI (see p. 51) of the Y array, 
or it is calculated from values supplied by the user of option 32 in P(9) 
and P(10) as FY = P(10) • 10 P ^ -6 . It is used to calculate KY and the 
values of YLA BEL. 
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FYLAB 


FYLAB5 

I 

n 

IK 

IL 

ILIM 

IM 

IMIN 

IXLAB5 

IYLAB5 

JX 

K 

KFD 


Variable FORMAT array with which the line of y grid-labels is written. 
It is of the form 


(2H , 20X, llFlO.d) 

where 0 < d < 6 . A value of d is inserted by the program (as KFD) 
after the call to PISTUG to process the Y array. See KFD, page 52. 

Initial value for FYLAB (5) . 

Index of many loops . 

Inner index of double DO loop that initializes the Y and K arrays for 
DUPY. 

Counter used as a subscript in the loop that associates the correct 
plotting character with the current value of IMIN for DUPY and NO DUP . 

Index of the loop that finds the correct plotting character for DUPY and 
NO DUP. 

Number of executions of the loop that finds XMIN and IMIN. 

Index of the loop that calculates KY, finds the associated plotting char- 
acter, and positions it in the A array. 

Value of I for the X(I) that is the smallest value in the X array at any 
particular time. 

Equivalent to XLAB5. 

Equivalent to FYLAB(5) . 

Initialized to 1, increased by 1 if option 64 is being used and a point has 
been placed in the line, increased by 2 if the line to be printed has an 
x-label, this switch value may be 1, 2, 3, or 4. It is used to control 
the selection of the correct WRITE statement for each line. 

Array name, the third argument in the call list of this subroutine. The 
first element is the sum of the option numbers being used; the second is 
the number of curves. The remaining elements of the array are pre- 
scribed by the Variation and the options being used. 

See COMMON/JOLO/ - Dictionary, page 52. 
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KLAST 

K 

KN 

KODE 

KOUT 

KOUTX 

KPC 

KPCSTD 

KPWR 

KQUIT 

KSW64 

KSWI 

KSX 

KSY 

KSYLAB 

KTIMES 


Limit of the loop that finds the correct plotting character for DUPY and 
NO DUP. 

Array name, the third in the call list of this subroutine. Its contents 
are KODE, KN, and other information as required by the Variation and 
options being used . 

Number of curves to be plotted. It is supplied by the programmer in 
K(2) . 

Sum of the option numbers. It is supplied by the programmer in K(l). 
Branching on single bit positions using MASK!, MASK2, MASK4, 

MASK8, MASK16, MASK32, and MASK64 serves to identify the options 
being used. Option 8 has been discontinued but a message is printed if 
MASK8 causes a branch. 

Switch initialized to 1 and set to 2 after all points have been processed.- 

Value in the special heading of option 64 for the x coordinates. It is the 
negative of KPWR of the X array . 

Equivalent to TPC. Temporary storage for a plotting character before 
it is positioned in the A array. 

Array containing the Six standard plotting characters (* + 0 X = □). 

See COMMON/JOLO/ - Dictionary, page 52. 

Switch initialized to 1. It is set to 2 just before the last line (a labeled 
x grid-line) of the plot is printed. 

See COMMON/JOUD/ - Dictionary, page 52. 

Switch initialized to 1 and set to 2 if Variation I (DUPX) is used. 

Fixed point form of the value supplied in P(6) by the programmer using 
option 16. P(7) and P(8) are multiplied by io^SX-8 to get FX and DX. 

Fixed point form of the value supplied in P(9) by the programmer using 
option 32. P(10) and P(ll) are multiplied by 10 KSY " 8 to get FY and DY. 

Switch initialized to 1 and set to 2 after the first line of y grid-labels 
has been written. 

Value of KN - 1. The number of executions of the loop that stuffs the Y 
and K arrays for DUPY. 
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KTL 

KX 

KY 

KYL 

L 

LA BOUT 
LCTR 

LL 

LS 

M 

MASK1 

MASK2 

MASK4 

MASK8 

MASK16 

MASK32 

MASK64 


Number of executions of the loop whose index is IM. It is 1 for DUPY 
and NO DUP and KN for DUPX. 

Line number of a particular point. It is computed from px. - FX)/DXj 
rounded. The first x grid-line is line number zero. 

"Printing position" of a particular point. It is computed from 
[(y^ - FY)/DY] rounded. The first y grid-line character in "printing 
position" zero is actually the second element of the A array. 

Value of KY + 2 if the y-coordinate falls on the plot; -1, if KY is nega- 
tive; 104, if KY > 101. 

Subscript computed and used in the loop that stuffs the Y and K arrays 
for DUPY. 

See COMMON/JOLO/ - Dictionary, page 52. 

Line counter. It is set to zero for the first x grid-line and is stepped by 
1 after each line is printed. 

Subscript computed and used in the loop that computes each KY and finds 
its associated plotting character. 

Logical variable set to TRUE whenever a point has been found for the 
current line. If option 64 is in effect, LS is used to control a branch 
(to modify JX) and is set to FALSE immediately after the branch. 

Frequency of the x grid-lines. It is either set equal to 10 or supplied in 
P(3) by the programmer if option 2 is being used. If the programmer 
specifies zero, M is set to 1000, and only the first and last x grid-lines 
will be printed. 


VSee KODE, page 41. 
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MM 

N 

NCTR 

NPTS 

NPTST 

NY 

P 

PWR10X 

PWRIOY 

RMARK 

STUG 

TDY 

TEMP 

TLINX 

TONLY 

TPC 

X 


Subscript computed and used in the loop that stuffs the Y and K arrays if 
DUPY is used. 

See COMMON/JOLO/ - Dictionary, page 53 . 

Counter of the number of x -values that have been processed. It is ini- 
tialized to 1 and tested for equality to ILIM. It is increased by 1 each 
time the test fails. 

Number of points in the X array (or Y array) supplied by the program- 
mer using DUPX (or DUPY) in K(3). 

Sum of the numbers of points in all the curves (supplied by the program- 
mer in K(3), K(5), . . . K(2 * KN + 1)) for NO DUP. It is KN multi- 
plied by NPTS for DUPY. 

Frequency of the y grid-lines. It is either set to 10 or supplied in P(4) 
by the programmer using option 4. If it is supplied as zero, it is re- 
placed by 1000 and only the left-most y grid-line will be printed. 

Array name, the fourth argument in the call list of this subroutine. Its 
contents are prescribed by the Variation and the options being used. 

Value of 10 KSX " 6 . 

Value of 10 KSY_6 . 


Out-of-range character (=) . 

See COMMON/JOLO/ - Dictionary, page 53. 

Intermediate storage for the value 10* DY, used in the calculation of the 
Y LA BEL array. 

Temporary storage. 

See COMMON/JOLO/ - Dictionary, page 54. 

See COMMON/JOLO/ - Dictionary, page 54. 

Current plotting character . Equivalent to KPC . 

Array name, the first argument in the call list of this subroutine, con- 
taining the values of the variable to be plotted on the vertical scale 
(down the page). The contents are destroyed during execution. 


PLOTMY 


43 



XGL 

XGRID 

XLAB 


XLAB5 

XLAB 51 
XLAB64 


XLABEL 

XMIN 

XOUT 

XYX 

Y 

YGRID 

YLABEL 


Logical variable set to TRUE whenever the value of LCTR is evenly 
divisible by M. 

x grid-line character - a minus sign. 

Variable FORMAT array used to print a line when JX = 3 (see p. 40). 
It is of the form 


(2H , 18X, F9.d, 104A1) 

where 0 < d < 6. The value of d is inserted by the program (using 
XLAB5) after the call to PISTUG to process the X array. See KFD, 
page 52. 

A word initialized to F9. 0 and altered to F9. d by using KFD (see p. 52) . 
XLAB5 is then moved into the FORMAT arrays (XLAB and XLAB64) 
used for lines that have an x-label. It is equivalenced to IXLAB5. 

Initial value for XLAB5. 


Variable FORMAT array used to print a line when JX = 4 (see p. 40). 
It is of the form 


(2H , 12X, F6.3, F9.d, 104A1) 

where 0 < d < 6. The value Of d is inserted by the program (using 
XLAB5) after the call to PISTUG to process the X array. See KFD, 
page 52. 

Current value of the x label. It is computed only for every tenth line 
from XLABEL = FX + LCTR-DX. 

Smallest value in the X array at any particular time. 

Value of modified by F10X for option 64 printouts. 

See COMMON/JOLO/ - Dictionary, page 54. 

Array name, the second argument in the call list of this subroutine, 
containing the values of the variable to be plotted on the horizontal scale 
(across the page) . 

y grid-line character, the digit 1 . 

Array holding the 11 labels for the y axis. 
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m. PISTUG 


A. Description 

This routine is internal to the plotting system and is used by PLOTXY and 
PLOT MY. The calling statement is CALL PISTUG (ARRAY). ARRAY is an array in 
any order . All other information is transmitted through 

COMMON/JOLO/F, DX, TUNX, N, LABOUT, KPWR, KFD 
COMMON/JOLO/XYX, FORY, STUG, TONLY, KSW64 

PISTUG uses the minimum and maximum values of the array to compute the total 
range. When scaling parameters are being computed for the y scale, 101 print posi- 
tions are available to cover the range. For the x scale, an arbitrarily chosen number 
of lines is assigned. In either case, the scale factor computed is always a value of 
DxlO n , where D is a member of the set [2, 2. 5, 5, 10] . 
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B. Program Listing 


SUBROUTINE P I STUG ( ARRAY ) 

LOGICAL XYX,FORY, STUG,TONLY ,XGL ,LS,KSW64 
DIMENSION ARRAY! 1 ) 

COMMON/ JOLO/F,DX,TL INX,N,L ABOUT, KPWP,KFD 
COMMON/ JOLO/X YX, FORY, STUG, TONLY, KSW64 
KHAR(XMAX) = I NT ( ALOGt XMAX ) /2. 302585+40 .0 ) -40 
126 XI = ARRAY (1 ) 

12 8 I FT X YX ) GO TO 133 
130 DO 132 J = 2 , N 

132 XI = AM INI (XI, ARRAY < J ) > 

133 IF (STUG) XI =F 

134 XN - 0.0 

136 DO 146 J = 1 , N 

138 DIF = ABS ( XI— ARRAY ( J ) ) 

140 IF( DIF.LE.XNIGO TO 146 
142 XN=D I F 
144 I HOL D=J 

146 CONTINUE 

147 XN = ARRAY (I HOLD) 

148 IF(KSW64 ) KPWR = KHAR ( AMAX1 ( ABS ( XI ) ,ABS ( XN) ) ) 

149 I F ( T ONLY ) GO TO 240 

150 T L I N= 1 01 . 

152 I F { .NOT. FORY) TL IN = TLINX 
154 C5 = (XN-XD/TLIN 
156 C6 = ABS ( C5 ) 

158 IF(C6.E0.0. ) GO TO 300 

159 K7 = KHAR ( C6 ) 

160 C8 = 1 0 .**K7 
162 C9 = C6/C8 

164 IF( (2.5-C9) .LE.O.O) GO TO 172 
166 0=2 . 

168 IF( (2.0-C9) .LE.O.) D=2.5 
170 GO TO 176 
172 D=5 . 

174 IF( (5. -C9). LE.O.O) D=10. 

176 Cll = 0*08 

178 DX = S I GN ( C 1 1 ,C 5 ) 

179 HUND = 100. *DX 
?40 K7 = KHAR ( ABS (DX) ) 

250 KFD = 0 

252 I F ( K7 ) 260,270,254 
254 IF(K7.GE.5) LAB0UT=2 
256 GO TO 270 
260 KFD = 6 

262 IF(K7.LT. (-7) ) LABOUT = 2 
264 IF(K7.GT.(-6) ) KFD = -K7 
270 I F ( STUG ) GO TO 230 
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I NT ( A B S ( X 1 ) / C 1 1 ) 


182 

184 

186 

188 

190 

192 

193 

194 

195 

196 

198 

199 

200 
202 
204 
206 
208 
210 
212 
214 
220 
222 
224 
226 
228 

229 

230 

300 


KC1 2 « 

JJ = 1 
IF(Xl) 188* 192 * 190 
JJ = 3 

I F ( DXVLT *0 • ) J J = JJ + 1 
GO TO (I93 t 194,195,196), JJ 
KC 1 4 = KC 12 
GO TO 198 
KC 14 * KC 12+1 
GO TO 198 
KC 1 4 = -KC12-1 
GO TO 198 
KC14 = -KC 12 
KC 1 3 = MOD (KC 12 ,10) 

KC 1 5 = KC12-KC13 
KC 1 8 = KC 1 5 

GOTO (212 ,202,202,210) , jj 
KC 1 8 = KC 1 8+1 0 
I F ( KC13 »NE #9 ) GO TO 210 
KC 1 8 = KC 1 4 
GO TO 212 

I F ( X 1 * L T * 0 • ) K C 1 8 = -KC 1 8 
F=C 1 1 AFLOAT (KC18) 

I F ( • NOT • FOR Y ) GO TO 230 
TEMP = F+HUND 

GO TO ( 22.4 » 22 8 1 224 1 22 8 ) * JJ 


IF(TEMP.G£.XN> 

GO 

TO 

230 

GO TO 229 




IF ( TEMP.LE.XN) 

GO 

TO 

230 

F=C1 1 AFLOAT (KC 

141 




CONTINUE 
RETURN 
DX = 0 • 

LAB OUT =3 
GO TO 230 
END 
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C . Block Diagram 

PISTUG 


150 

[TUN = 10l|\ 


XI = A(l), 
AMIN, or F 
XN=A(IH0LD) 


Is scale / 
factor (DX)/ No 


126 

"I XI- A(l) | 


Scale x / 
for / 

\plqtxy/ No 


Yes 

VI 

XI = AMIN 

133 , 




Is 16 (or 32) y 
Vbeing used/ 


Calculate maximum 
possible scale factor 
= C5 and |C5| 


1 59-162 V™ 

I K7, C8 C9 1 


Set signal for 
error message 
DX - 0 
LA BOUT = 3 


\150-179| 

Calculate 

< Error HUND and 

exit scale factor 

DX 


( Find index 
(IHOLD) 
of element 
furthest 
L from XI 

147 I 
k XN= A(IHOLD) 


^176-179 / 

)X|. DX, HUND 




Test for acceptable 
label size and find 
KFD ( M d M for for- 
mat specification) 


Is 16 (or 32)/ V 
\being used / ^ es \ 


Compute setting of 
switch JJ and best 
possible value of F 
(first label) 


250 | 

I Set KFD = "d" = 0 


260 


Set KFD 

. "d" = & 

262 i 



K7 too / 

V I'ttle / F~~. “ 

\ ? / Yes Se signa or 

\ / ► error message 

Y (error) labOUT - 2 


ig / Use KFD = 0 
7 / Yes and set signal 
/error)* for error mes- 
/ sage LA BOUT =2 

No 


220 

I Calculate last label 
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p. PISTUG Dictionary 

C5 Value of the maximum possible scale-factor. It is found by dividing 

the range XN - XI by TUN. 

C6 Absolute value of C5. It may not be zero. If it is, the error switch 

LABOUT is set to 3 and the plot is terminated. It is used to deter- 
mine C9. 

C8 Intermediate storage used in the calculation of C9 and DX. 

C9 When the maximum possible scale-factor C5 is expressed as x.xxxx 

multiplied by 10 n , C9 is the absolute value of the x.xxxx factor. It is 
used to select D, the next larger value from the set of permissible 
values [2, 2.5, 5, lOj. 

Cll Intermediate storage used in the calculation of DX and of F. 

D Smallest member of the set [2, 2. 5, 5, 10] that is larger than C9. It 

is used when PISTUG is calculating the scaling parameter DX. 

DIF Absolute value of the distance between XI and the element furthest 

from XI. 

DX See COMMON/JOLO/ - Dictionary, page 52. 

F See COMMON/JOLO/ - Dictionary, page 52. 

FORY See COMMON/JOLO/ - Dictionary, page 52. 

HUND One hundred times the scaling parameter DX. 

IHOLD Index of the element in the array which is furthest from XI. 

J Index of the search loop in which DIF and IHOLD are found. 

JJ Switch used to control the calculation of the starting value, F, 

according to the algebraic signs of XI and DX. 

K7 Characteristic of the log^Q of the absolute value of the scale -factor 

DX. It controls the value of KFD, the number of decimal places in 
the grid-label, as follows: 
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KC12 

KC13 

KC14 

KC15 

KC18 

KFD 

KHAR(Z) 

KPWR 
KSW64 
LA BOUT 
N 

STUG 

TEMP 

TUN 

TUNX 

TONLY 


KFD = 0 if K7 is + or zero. 

KFD = -K7 if K7 is negative and greater than -6. 

KFD = +8 if K7 is negative and equal to or less than -6. 

If K7 is greater than 4 or less than -7, the error switch LABOUT is 
set to 2 to cause an error message to be printed below the plot. 

The symbol K7 was inadvertently also used for intermediate storage of 
the characteristic of log 10 C5 during the calculation of C9. 

"\ 


Intermediate storage used in the calculation of the "best possible" 
I value of F. 


See COMMON/JOLO/ - Dictionary, page 52. 


Arithmetic function that computes the integer characteristic of 
loglO Z- 

See COMMON/JOLO/ - Dictionary, page 52. 

See COMMON/JOLO/ - Dictionary, page 52. 

See COMMON/JOLO/ - Dictionary, page 52. 

See COMMON/JOLO/ - Dictionary, page 53. 

See COMMON/JOLO/ - Dictionary, page 53. 

Temporary storage. 

When PISTUG is scaling a Y array, TUN is set equal to 101. When 
an X array is being scaled, TUN is supplied by the calling routine as 
TUNX. TUN is used to compute C5. 

See COMMON/JOLO/ - Dictionary, page 54. 

See COMMON/JOLO/ - Dictionary, page 54. 
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XI Minimum or maximum value of the array being scaled. It is set equal 

to the first element of the array when PISTUG is scaling an X array 
for PLOTXY, otherwise it contains the minimum value of the array 
after a search. If STUG is TRUE, XI is set equal to the scaling- 
parameter F supplied by the calling routine. XI is used to compute 
C5 and KPWR, and several times in the calculation of F. 

XMAX Dummy argument of arithmetic function KHAR, page 50. 

XN Minimum or maximum value of the array being scaled. It is set equal 

to zero and replaced by the element furthest away from XI. It is used 
to compute C5 and in the calculation of KPWR. 

XYX See COMMQN/JOLO/ - Dictionary, page 54. 
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E . COMMON/JO LO/ - Dictionary 


DX 

F 

FORY 

KFD 

KPWR 

KSW64 

LABOUT 

1 

2 


Scaling parameter for the scale factor for either scale . It may be cal- 
culated by PISTUG or supplied indirectly by the user. For PLOTMY, 
DX for the x scale must be positive. 


Scaling parameter for the starting-value for either scale. It may be 
calculated by PISTUG or supplied indirectly by the user. 


Logical switch set to TRUE by PLOTXY or PLOTMY only if the call 
to PISTUG is for the y scale. 


Number of decimal places required in the grid-label for the array be- 
ing processed by PISTUG. It is used to modify the variable FORMATS 
(FXLABS, FXLABM, FYLAB in PLOTXY and XLAB, XLAB64, FYLAB 
in PLOTMY) used to write the grid-labels. This modification is made 


by adding KFD to a word initialized to blank, blank, F9. 0 


Characteristic of the log^Q of the maximum of the absolute values of 
XI and XN. (The minimum and maximum of the array being pro- 
cessed by PISTUG. ) It is calculated by PISTUG only if option 64 is 
being used. It is used by PLOTXY and PLOTMY to prepare the coor- 
dinates being printed at the left of the plot so they will fit a F6. 3 
FORMAT specification. A special heading is written displaying the 
conversion factors. 


Logical variable set to FALSE by PLOTXY and PLOTMY and set to 
TRUE if option 64 is used. It is used by PISTUG to control calculation 
of KPWR and in PLOTXY and PLOTMY to control printing of coor- 
dinates at the left of the plot. 

This is a computed-go-to index used to control branching to normal 
termination or to various error terminations. 

Normal ending of a plot; a blank record. 

Error message BAD LABELS is written below the plot. This is an 
indication that one or more x or y grid-labels are either too large 
or too small to be printed correctly in the assigned label field 
(F9.d - where 0 < d < 6). This is checked by PISTUG for F and by 
PLOTXY and PLOTMY for successive labels. The relative positions 
of the plotted points are always correct; the labels may or may not be. 
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3 The plot is terminated immediately if: 

(a) The range of an array being scaled by PISTUG is found to be zero. 
The error output is 

PLOTXY N. G. followed by X(l), Y(l), X(2), Y(2), K(l), P(l) 

PLOT MY N.G. followed by X(l), Y(l), X(2), Y(2), K(l), K(2), K(3) 

(b) The values in the X array are not in monotonic order when 
PLOTXY is used. The error output is X OUT OF ORDER followed 
by the message in (a) . 

4 The plot is terminated immediately if option 16 is used in PLOTXY or 
PLOTMY and a value of x is found outside the user-specified starting- 
value (Rewind Carriage Error) . The error output is the same as in (b) . 

5 If the search in PLOTMY for the plotting character fails, the plot is 
terminated and the error message ERROR IN K ARRAY is written. 

N Number of points in the array sent to PISTUG to be scaled. If the call 

is from PLOTXY, N has been supplied by the user in P(l). If the call 
is from PLOTMY: 

Variation I - N is the number of points in the X array (supplied 
by the user in K(3)). 

Variation II - N is the number of points in the X array. It is cal- 
culated by multiplying the number of points in the Y array 
(supplied by the user in K(3)) by the number of curves, KN 
(supplied by the user in K(2)) . 

Variation III - N is the total number of points to be plotted and is 
calculated by adding the number of points for each curve 
(supplied by the user in K(3), K(5), K(7), etc). 

STUG Logical switch set to TRUE by PLOTXY or PLOTMY when PISTUG is 
not required to calculate F. 
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TLINX 

TONLY 

XYX 


If PISTUG is scaling an X array, the calling routine must supply the 
number of lines into which the plot must be fitted. The empirically 


chosen function of N, TLINX = 
35 points onto one page. 



fits most plots of less than 


Logical switch set to TRUE by PLOTXY or PLOTMY if PISTUG is not 
required to calculate F or D. 

Logical switch set to TRUE only if the call to PISTUG is for the x 
scale from PLOTXY. 
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IV. SORTXY 


A. Description 


CALL SORTXY (V, W, NPTS) 

This subroutine rearranges the NPTS elements of the V array in order of in- 
creasing size. The elements of the W array are moved to maintain the original pair- 
relation; that is, if the fifth element of the V array is moved to the first position of 
V, the fifth element of W is moved to the first position of W. 

B. Program Listing 


SUBROUTINE SORT XY ( X , Y ,NPTS ) 
DIMENSION X ( 1 ) , Y ( 1 ) 

100 N=NPTS 
102 NN=N-1 
104 DO 140 KT=1,NN 
XMI N=X (KT ) 

J AD=KT 
JKI_ = KT + 1 

112 DO 120 JK=JKl,N 

114 IF ( XMIN-X ( JK) ) 120,120,116 

116 XMIN=X(JK) 

118 J A D - J K 
120 CONTINUE 
122 YMI N=Y ( JAD ) 

X ( J AD) = X ( KT ) 

Y ( JAD ) = Y ( KT ) 

X ( KT ) = XMIN 
Y ( KT ) = YMI N 
140 CONTINUE 
RETURN 
END 
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V. SKALE 


A. Description 


CALL SKALE (NPTS, A, KRSTR) 

This subroutine finds the largest absolute value of the NPTS elements of A, and 
computes the characteristic of its log^. If the characteristic K is -2 < K < 4, 
KRSTR is set to zero and control returns to the calling program. If K > 4 or <-2, 
each element of A is multiplied by a power (KRSTR) of 10 to transform the array to 
suit PLOTXY and PLOTMY. 

KRSTR is returned to enable the user to ReSToRe the array or record how it has 
been altered, 

B. Program Listing 


SUBROUTINE SKALE ( NPTS, X, KRSTR ) 

DIMENSION X ( 1 ) 

100 N-NPTS 
102 XMAX=X { 1 ) 

104 DO 106 J = 1 1 N 

106 XMAX = AMAX1 (ABS(XMAX) , ABS I X ( J ) ) ) 

108 KHAR = I NT ( AL0S ( XMAX ) /2. 302585+40 »0 )-40 
116 I F ( ( 4-KHAR ) * ( KHAR+2 ) ) 120 ».l 1 8 » 118 

118 K PWR=0 

119 GO TO 130 

120 KPWR=3-KH AR 

122 FACT=10.0**KPWR 
126 DO 128 J=1,N 
128 X(J)=X(J)*FACT 
130 K R S T R - K P W R 
140 RETURN 
END 
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VI. AND 


A. Description 


FORTRAN G compilers do not provide the logical AND function available in 
FORTRAN H. The following FORTRAN function subprogram will satisfy the require- 
ments of PLOTXY and PLOTMY. Despite its name, it is not usable as the AND func- 
tion for other purposes. 


B. Program Listing 


FUNCTION AND (KODE,MASK) 

100 K = KODE 

IF(KODE.EQ.O) GO TO 110 

IF ( MOD (K/ MASK, 2) .EQ.O) GO TO 110 

AND= 1. 

RETURN 
110 AND= 0. 

RETURN 

END 


57 



VII. INPUT -OUTPUT CONSIDERATIONS 




The printer plotting system used at Lewis operates under OS with no print con- 
trol. All records (limited to 132 characters) written with a WRITE(6,###) state- 
ment are printed immediately . If a supervisory program is to control printing, 
records written by the plotting routines must be treated somewhat differently than 
ordinary output records, for example, automatic page-spacing must be disabled. For 
this reason, all FORMAT statements begin with a field of two blanks. This is pro- 
vided to permit easy modification of the FORMAT statements if the supervisory pro- 
gram requires identification within each record. 

For convenience in making changes the WRITE and FORMAT statement numbers 
for PLOTXY and PLOT MY are listed below. 


WRITE 

FORMAT 


100 

500 

Writes a blank record. This causes a blank 
line immediately preceding the plot. 

260 

502 

Writes a special heading to the left of the 
plot if option 64 is used. 

302 

FYLAB 

This FORMAT statement is used to write 
the lines of y grid-labels and is of the gen- 
eral form (2H , 20X, llFlO.d). It is an 

array to permit setting of the value of "d" 
(in the F specification) to suit the range of 
the Y array. See KFD, page 52. 

702 

506 


704 

520 

>• These are error messages. 

706 

508 

J 


710 

510 

This FORMAT statement is used to write a 


blank record following the second line of y 
grid-labels. 


712 512 This is an error message if the search for 

(PLOTMY ONLY) the associated plotting character fails (see 

D, top of p. 14). 
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WRITE 

6250 

6260 

6270 


6280 


FORMAT 

504 

505 

PLOTXY FXLABS 
PLOT MY XLAB 


This FORMAT statement is used to print a 
line when it has no x -label nor option 64 out- 
put to the left of the plot. 

This FORMAT statement is used to print a 
line that has no x-label but does have option 
64 output to the left of the plot. 

This FORMAT statement is used to print a 
line containing an x-label but no option 64 
output to the left of the plot. It is of the gen- 
eral form 

(2H , 18X, F9 . d, 104A1) 

Using an array permits setting the value of 
”d M (in the F specification) to suit the range 
of the X array. See KFD, page 52. 


PLOTXY FXLABM This FORMAT statement is used to print a 

PLOT MY XLAB64 line containing both an x-label and option 64 

output to the left of the plot. It is of the gen- 
eral form 


PLOTXY (2H , 4X, 2F7.3, F9.d, 104A1) 

PLOT MY (2H , 12X, F6.3, F9.d, 104A1) 

Using an array permits setting the value of 
"d" (in the F specification) to suit the range 
of the X array. See KFD, page 52. 
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SAMPLE PLOTS 
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elements in the array that is the first argument (2) KODE a 69 a (64 + 4 + 1), the desired plotting character in P(l), and P(3) * 20. 
The vertical gridline frequency is changed, the coordinates are printed to the left of the plot and the plotting character is the pro- 
grammer 1 s choice. 



*EXP(-3.62-Y*YM.<> 
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FOUR VALUES OF THE DE PE NOE NT VARIABLE ( ACROSS ) 
FOR EACH VALUE OF THE INDEPENDENT VARIABLE (DOWN 
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Figure 3. - This plot displays the use of PLOTMY - Variation I, KODE 



FOUR VALUES OF THE DEPENDENT VARIABLE (DOWN) 

FOR EACH VALUE OF THE INDEPENDENT VARIABLE IACROSS1 
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Figure 4. - This plot displays the use of PLOTMY - Variation II, KODE = 0. The same data are plotted here as in figure 3; 
only the Variation number and the order of the first two arguments have been changed. 
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"T&? aeronautical and space activities of the United States shall be 
conducted so as to contribute , . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space , The Administration 
shall provide for the widest ptdgticable and appropriate dissemination 
of information concerning its activities and the results thereof/’ 

— National Aeronautics and Space Act of 1958 
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NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


' 


TECHNICAL REPORTS: Scientific and TECHNICAL TRANSLATIONS: Information 

|H£ technical information considered important, published in a foreign language considered 

complete, and a lasting contribution to existing to merit NASA distribution in English, 

jl / knowledge. ' M 




TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: : ; 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 
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Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 
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